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The chemistry of planetary nebulae in the outer regions of M31 
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ABSTRACT 

We present spectroscopy of nine planetary nebulae (PNe) in the outskirts of M31, all but 
one obtained with the 10.4 m GTC telescope. These sources extend our previous study of the 
oxygen abundance gradient of M31 to galactocentric radii as large as 100 kpc. None of the targets 
are bona fide members of a classical, metal-poor and ancient halo. Two of the outermost PNe 
have solar oxygen abundances, as well as radial velocities consistent with the kinematics of the 
extended disk of M31. The other PNe have a slightly lower oxygen content ([0/H]~—0.4) and 
in some cases large deviations from the disk kinematics. These PNe support the current view 
that the external regions of M31 are the result of a complex interaction and merger process, 
with evidence for a widespread population of solar-metallicity stars produced in a starburst that 
occurred ~2 Gyr ago. 

Subject headings: galaxies: abundances - galaxies: individual (M31) - ISM: abundances - planetary 
nebulae: general 


1. Introduction 


Planetary nebulae (PNe) are valuable tracers of 
stellar populations in all types of galaxies. They 
are widely used as distance indicators, via the in¬ 
variant bright cutoff of the PN lum inosity function 
(PNLF; e.g. ICiardullo et all 2002h, an d as tracers 
of the luminosity (Buzzoni. Arnaboldi, & Corradi 


l2006h . dynamic s (e.g. Merrett et al. 200 611 and 


chem istry (e.g. iKwitter et al. l2012t iBalick et al 


120131 . hereafter papers I and II, respectively) of 
their stellar progenitors. 

Recently, PNe have been used to determine the 
metallicity gradient_in the disk o f the nearby sp i¬ 
ral galaxies M33 ( Magrini. Stanghellini fc Villaven 
120091) and M31 (papers I and II). These gradi¬ 
ents provide direct constraints to models of disk 
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formation and evolution. In particular, measure¬ 
ments in the outer regions of the disks have the 
power to test novel ideas such as the importance of 
external distu rbances to p r oduc e extended young 
stellar disks ( Werk et al.l l201ll ). In the case of 


M31, it has been proposed that a tidal encounter 
with M33 that occurred about 3 Gyr ago would 
have produced a vas t extended disk with homo¬ 
geneous metallicity ( Bernard et al.1 20121 2015h . 
Such an extended disk has been reported before 
( Ibata et~ahl 12005b and PNe provide the oppor¬ 
tunity to test its predicted chemical content. In 
earlier papers, we have studied PNe out to depro- 
jected galactocentric distances of 60 kpc, finding 
nearly solar O/H abundances and a flat O/H gra¬ 
dient in these regions. In this work, we extend our 
study by obtaining high-quality spectra of another 
nine PNe at larger radii in the outskirts of M31. 


2. Observations 


The target PNe are listed in Table [l] All but 


one were selected from the list of iMerrett et al 


( 2006h to be at large galactocentric distances and 
sufficiently bright to obtain high quality spectra 
with a 10m teles cope. The putat i ve hal o PN M31- 
372 observed by Ijacobv fc Ford ( 19861) was addi- 
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Table 1 

Basic properties of the target PNe 


Name 

R.A. a 

[J2000] 

Dec a 

[J2000] 

m(5007) a 

[mag] 

£ 

[deg] 

V 

[deg] 

dapp 

[deg] 

dapp 

[kpc] 

ddisk 

[kpc] 

V©, sys a 

[km s -1 ] 

V*// 

[km s -1 ] 

M2507 

00 48 27.2 

+39 55 34.3 

21.23 

1.10 

-1.33 

1.73 

23.2 

106.3 

-147 

180 

M2538 

00 36 28.8 

+39 35 26.4 

20.25 

-1.21 

-1.67 

2.06 

27.7 

28.0 

-426 

no 

M2539 

00 36 12.6 

+39 35 41.9 

21.16 

-1.26 

-1.66 

2.08 

28.0 

28.0 

-426 

no 

M2541 

00 35 09.1 

+39 28 25.2 

21.78 

-1.46 

-1.78 

2.30 

31.0 

31.3 

-456 

80 

M2543 

00 35 50.7 

+42 21 04.5 

21.61 

-1.27 

1.09 

1.68 

22.6 

105.8 

-272 

40 

M2549 

00 36 27.2 

+42 06 21.9 

21.35 

-1.17 

0.85 

1.44 

19.4 

90.9 

-17 

290 

M2566 

00 49 28.3 

+40 59 53.9 

21.26 

1.27 

-0.26 

1.30 

17.4 

73.8 

-247 

10 

M2988 

00 52 00.0 

+43 03 23.5 

21.98 

1.69 

1.81 

2.48 

33.3 

36.2 

-98 

20 

M31-372 b 

00 46 41.5 

+43 59 03.7 

22.6 

0.71 

2.72 

2.81 

37.8 

77.6 

-60 

100 


a From IMerrett et alj ll2006l '). except for M31-372. 

b Coordinates are from fjacobv fc Fordl J I98f| ). V sys (uncertainty 40 km s —1 ) and m(5007) (uncertainty 0.3 mag) 
are estimated from our spectrum. 

Note. —The following p arameters for M31 have been assumed throughout the paper: distance 770 kp c 
dFreedman fc Madord[l99of); center at RA(J20 00)= 0 42 44.3 and Dec(J2000 )=-h41 16 09.0 (IMerrett et £>11120061) ; 
disk inclination i=77 .7 JPe Va ucouleur^[l956 ) and position angle PA=37°.7 (IMerrett et al. IM); heliocentric sys¬ 
temic velocity Van.=—309 kms llMerrettet*lrL l200fiT l. 


tionally included. 

Spectroscopy of eight of these PNe was obtained 
in different nights in October 2013 in service-queue 
mode at the 10.4m GTC telescope on the island 
of La Palma, Spain. The OSIRIS instrument was 
used in its longslit mode. The combination of 
grism R1000B and a slit width of 0".8 provides 
a spectral dispersion of 0.21 nm per (binned x2) 
pixel, a resolution of 0.63 nm, and a spectral cov¬ 
erage from 370 to 785 nm. Seeing varied between 
0".6 and l".l (full width at half maximum), and 
spatial scale along the slit was 0".254 per binned 
pixel. Data were generally obtained under pho¬ 
tometric weather conditions and grey moon. The 
slit was oriented along the parallactic angle for 
all targets and standard stars. Total exposure 
times per target varied from 120 to 135 min, de¬ 
pending on the brightness of the source, split into 
three or four sub-exposures. The GTC calibration 
plan provides, in each observing night, at least one 
spectrophotometric standard to be used for flux 
calibration. 

Additionally, the PN M2538 was observed with 
the Dual-Imaging Spectrograph at Apache Point 
Observatory (New Mexico, USA) in 2011 Octo¬ 
ber. We observed through a 2" slit oriented along 


the parallactic angle. The B1000 grating provided 
coverage from 370 to 505 nm at 0.25 nm reso¬ 
lution; the R300 grating covered 520-960 nm at 
0.70 nm resolution. Standard bias, flatfield, and 
emission-lamp exposures were taken and applied 
to the data. Flux calibration was obtained via ob¬ 
servations of the standard star BD+28 4211. We 
obtained six 20 min integrations of M2538, which 
were co-added after calibration. 

As PNe at the distance of M31 (770 kpc, see Ta¬ 
ble [D are spatially unresolved, 1 D spectra were 
extracted and reduced using twodspec in IRAF0. 


3. Physics and chemistry of the nebulae 


The spectra of the target PNe are displayed 
in Figure [l] Emission-line fluxes were measured 
by multi-Gaussian fit using splot. These fluxes 
are the input for the abundance determinations, 
carried out with EL SA, our five-level atom code 
( Johnson et al. 2006h . We used the same method 
as in Papers I and II, in order to produce a homo¬ 
geneous set of data and abundance estimations. 


1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni¬ 
versities for Research in Astronomy (AURA) under coop¬ 
erative agreement with the National Science Foundation. 
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Fig. 1.— The spectra of the target PNe. For the lower S/N spectrum of M2538, regions with high noise 
have been masked, and the whole spectrum was smoothed with a boxcar of three pixels. 









The reader is referr ed to those articles and to 
Milingo et al. ( 2010 1 for the details of the anal¬ 
ysis. The procedure includes corrections of the 
observed line fluxes for in terst ella r redd ening, us¬ 
ing the law of ISavaee fe Mathis! ( 1979h . and for 
the contamination of the hydrogen Balmer lines 
by coincident recombination lines of He ++ . Dered- 
dened fluxes were used to determine the electron 
temperature and density N e and T e with the ap¬ 
propriate line diagnostics. Note that T e {[O III]) 
was always properly determined, as the auroral 
line [O III] 4363 was accurately detected in all PNe 
(FigureQ]). The same applies for N e ([S II]), except 
for M31-372 and M2538 where the [S Il]6717,6731 
doublet could not be measured and a default N e of 
10000 cm -3 was assumed. In PNe where no direct 
low-ionization temperature diagnostic has been 
observed, the determination of the temperature to 
be used for calculating low-ionization abundances 
(T e ([N II])) depends on whether Heii4686 is ob¬ 
served. If it is, we ado pt the carefully derived 
result from llvaleil ( 1986ll that applies under this 
condition, i.e. T e ([N II] )=10300 K. If not, then 
we derive T e ([N II]) f rom T P (\O III]) acc ording to 


(1992). 


the prescription from lPagel et al. 

The detected lines allowed the calculation of the 
abundances of several O, Ar, N, S, and Ne ions. 
Total abundances were then determined using ion¬ 
ization cor recti on factors ( i cf) ca lculated as de¬ 
scribed in Kwitter fc Henry (200 1). The observed 
emission-line measurements, logarithmic redden¬ 
ing parameter c(H/3), T e and N e , ionic and to¬ 
tal abundances are listed in the tables in Ap¬ 
pendix A. They indicate that none of the PNe pre- 
sented_in this paper fulfill the o rigina l definition 
by Peimbert fe Torres-Peimbertl ( 19831 ) for Galac¬ 
tic Type I PNe of having either log (N/O) > —0.3 
or He/H> 0.125. This also applies to most of 
PNe studied in papers I and II, indicating at most 
a very moderate N and He enrichment for these 
bright PNe of M31. The PNe studied in this arti¬ 
cle also follow the relationships between O/H and 
other elemental abundances (such as N/O, Ne/H, 
and Ar/H) discussed in paper I. 


4. Location of targets within M31 

Figure [2] shows the position in the sky of 
the target PNe, including those discussed in pa¬ 
pers I and II. The geometric transformations 



|[deg] 

Fig. 2.— Location of the target PNe, overlaid on 
the DSS optical image. Circles and triangles in¬ 
dicate PNe with deviations from the model disk 
kinematics smaller and larger than 100 km s , 
respectively (see text). Red symbols are the PNe 
presented in this work, and black symbols those 
discussed in papers I and I I. Squares indicate th e 
HST/ACS fields studied bv lBernard et al.1 ( 2015 ). 
The ellipse indicates the R 25 radius of M31. 


of Huchra, Brodie & Kent (]199ll ). with the pa¬ 
rameters indicated in the Note of Tabic [T) were 
adopted to determine the R.A. and Dec offsets 
£ and 77 from the center of M31. These offsets, 
and the total apparent distance d app in the plane 
of the sky, are also listed in Table [I] The el¬ 
lipse in the figure marks the R 25 radius of M31 
(~30 kpc, see paper II). It corresponds to ~5 
scale-leng ths of the exp onential bright disk of M31 
( Walterbos fe Kennicuttl fl988l) . showing that the 
newly observed PNe are located well outside the 
familiar bright disk of M31. 

These outer regions of M31 are surprisingly 
complex. A huge halo with a radius >300 kpc 
and a smo oth star density d i stribu tion has been 


identified li lbata et al. 2007 . 2014I ). A number 


of structures, such as streams, loops, and over¬ 
density regio ns are seen in projection through - 
out the halo ( Lewis et al. 2013 ; Ibata et al. 2014 ). 
They are mostly associated with accretion of 
satellite dwarf galaxies, as expected in the stan- 
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dard hierarchical halo formation scenario. The 
most prominent structure is the so-called Giant 
Stellar Stream (GSS), which is likely the lat¬ 
est, most metal enhanced accretion event asso¬ 
ciated with the streams. In addition, outside the 
usual bright disk of M31, a smooth “exodisk” was 
found to extend out to a galactocentric distance 
of ~40 kpc, with detections as far as ~70 kpc 
( Ibata et al. 2009 1. A set of globular clusters in 
this zone forms a coherent kinematic system that 
i s similar to an extrap olation of M31’s inner disk 
( Velianoski et al.ll2014j h 

As far as metallicity is concerned, there is a 
significant overlap between the stellar [Fe/H] con¬ 
tent of the different structures in the outer re¬ 
gions of M31. Overall, the mean halo metal¬ 
licity decreases with radius from [Fe/H]~—0.7 
at R <30 kpc to [Fe/H]~—1.5 at R ~150 kpc, 
but the smooth halo component seems 0.2 dex 
more metal poor at all radii, and substructures 
span a wide range in metalli city covering r o ughly 
two orders of m agnitude (I Tanaka et al.1 l2010t 
Ibata et al. 2014th The strongest evidence of 


metal-rich stars in the halo (beyond the bound 
satell ite galaxies) at (|F e/Hl)~—0.5 is found in the 
GSS ( Ibata et al. 2014h . As with the extended ex¬ 
odisk, iBernard et al.l (l2015h find an average value 
of [Fe/H l~—0.3 for the several disk-like fields con¬ 
sidered. iGhaoman et al.l ((2006) find several fields 
in M31’s exodisk in which AGB stars ha ve rou ghly 
solar metal abundances. Figure 8 in ( Ibata et~ahl 


201411 illustrates the complexity of the metallicity 


distribution in the outer regions of M31 caused by 
the overlapping halo (with its ancient and more 
recent components) and the extended disk. 

Adopting the inclination and position angle of 
M31’s inner disk, the deprojected distances in the 
plane of the disk for the target PNe, ddisk, are 
indicated in Table |T| The radial velocities of 
PNe provide useful insights into their possible re- 
latio n to M31’s extended disk. They are adopted 
from Merrett et al.l ( 2006h . except for M31-372 for 
which a rough estimate could be obtained from 
our spectra (Table CO. We compared them with 
the average velocity expected at each position ac¬ 
cording to the kinem atic model of the extended 
disk presented bv llbata et al.1 ( 2005 ). The most 
significant differences (V*// >100 km s -1 , see 
Table [T]) are found in PNe M2507, M2538, M2539, 
M2549, and M31-372, suggesting that they may 
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Fig. 3.— The O/H abundance radial gradient. 
Symbols are like in Figure [2] Dashed lines are 
described in the text. 



instead be associated with the halo or some of its 
substructures. 

[Bernard et al. ( 2015h present colour-magnitude 
diagrams in fourteen HST/ACS fields sampling 
galactocentric distances similar to the PNe pre¬ 
sented in this work. They cover regions of the 
extended disk of M31, or some of the substruc¬ 
tures in the halo, none representing the metal- 
poor halo component. The positions in the sky 
of the HST fields are indicated by squares in Fig¬ 
ure [2] While some of these fields are not far from 
our PNe, no noteworthy association can be iden¬ 
tified, except for the case of M2538, M2539, and 
M2541 which seem related to the so-called Warp 
and G1 Clump, which are both disk components 
(but note the abovementioned deviations from the 
disk kinematics of two of these PNe). 


5. M31 O/H abundance gradient 

The O/H abundance radial gradient, compiled 
with all PNe studied in papers I and II (black 
symbols) and this work (red symbols), is pre¬ 
sented in Figure [3] It extends the metallicity gra¬ 
dient from this type of object to galactocentric 
distances larger than 100 kpc, assuming that all 
PNe are located in the plane of the disk of M31. 
Sources with radial velocities consistent with the 
disk kinematics (Vdi//<100 km s -1 , circles) are 
separated from those with lar ge deviation s from 
the disk kinematical model of Ibata et all ( 20051 
triangles). The graph shows that even at large 
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distances PNe with solar metallicity and disk-like 
kinematics exist. A least-square-fit to sources 
with Vdi//<100 km s -1 provides an almost flat 
metallicity slope ([0/H]=—0.0909—0.0004xd^isfc, 
with ([0/H])=—0.11). It is indicated by the up¬ 
per dashed line in the figure. As a note of cau¬ 
tion, the two PNe with solar oxygen abundance at 
ddisk > 60 kpc are M2543 and M2566. Both are 
located along the minor axis of M31’s disk. There¬ 
fore their deprojection factors to the disk plane 
are large, and radial velocities are close to the sys¬ 
temic velocity of the galaxy. They do not appear 
to be associated with any major sub structure in 
the halo (see also lMerrett et alJi2006 l. although it 
should also be noted that that the GSS pervades 
a large portion of the outer regions of M31 out 
to 100 kpc, and its debris extend beyond_the SE 


quadrant where it is more prominent (jlbata et al 
20141. 


Figure [3] also shows that the three PNe at 
ddisk > 60 kpc which deviate from the disk kine¬ 
matics have a lower metallicity (([0/H]}=— 0.37, 
lower dashed line). This indicates that they might 
be related to an older population, but their only 
slightly subsolar O/H seemingly precludes their 
association with the stars in the most metal-poor 
and anci e nt co m ponen t of the halo identified by 
e.g. Ibata et al. ( 2014 1. 

Before moving to a more general interpretation, 
two things should be noted. In this paper, we as¬ 
sume that the O/H abundance ratio in the PNe 
reflects the original oxygen content of the stel¬ 
lar progenitors, i.e. of the ISM from which these 
stars formed. There are some indication that oxy¬ 
gen may be enhanced in some AGB stars dur¬ 
ing the third dredge-up episodes, but this possi¬ 
ble effect is neglected h ere as no clear conclusion 
has been drawn yet ( Karakas fe Lattanzio| 201 it 
Henrv et al.ll2() 1 2t iDelgado-lnglada et al. 201 fil l. 


Also, it should be recalled that our target PNe 
were selected to be bright enough to allow an 
accurate chemical analysis. Indeed, all lie within 
2.4 ma g from the bright cutoff of the PNL F. Richer 


( 1993 1 in the LMC, and Jacoby & Giardnllo 


19991 1 in M31, suggested that there might be 


a slight dependence of metallicity with the [O III] 
luminosity, but concluded that the luminosities 
of the brightest PNe are nearly indepen d ent o f 
O/H. On the other hand, Magrini et al. ( 2004l l 
did not find any trend of chemical abundances 
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Fig. 4. 

The O/H abundances as a function of 


the absolute [O III] magnitudes of PNe in M31 
(top two panels, adopted distance 770 kpc) and in 
M33 (bottom panel, adopted distance 840 kpc). 
The dashed lines indicate the adopted extinction 
values (see text). For M31, symbols are like in 
Figure [2] 


as a function of the [O III] luminosity for PNe 
in the LMC and our Galaxy, concluding that the 
chemical abundances derived from the brightest 
PNe are representative of the total PN popula- 
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tion. Figure |4j compiled with our data of M31, 
shows some tendency of decreasing metallicity 
with [O III] luminosity. The upper panel displays 
our measured O/H abundances as a function of 
the absolute [O III] magnitude corrected for fore¬ 
ground and internal extinction using the c(H/3) 
values determined in our spectra. PNe in the 
brightest magnitude bin (M(5007)<—3.7) have 
some spread in metallicity, but on average have 
higher O/H abundances than PNe in the next 
magnitude bin. To minimize the effects of the - 
albeit shallow - abundance gradient, the middle 
panel confirms that the trend persists if sources 
are selected in a more limited range of galacto- 
centric distances. A similar relation is also seen 
seen in M33 flower panel), where we have used 
the data from Magrini, S tanghellini fc Villaver 
(2009), which were dereddened adopting the M33 
foreground and interna l extinction from Cepheids 


( Freedman et al. 200lh . The decrease of mean 


metallicity with the [O III] luminosity may ex¬ 
plain why the dispersion in the radial gradient of 
M31 at ddisk < 40 kpc is increased, compared to 
the corresponding graphs in papers I and II, with 
the addition of the slightly fainter PNe presented 
in this work. 

In a standard scenario where metallicity of a 
galaxy increases with time and the PN progeni¬ 
tors do not form or destroy oxygen, the observed 
correlation between luminosity and O/H con¬ 
tent is qualitatively consistent with the hypoth¬ 
esis that the most luminous PNe are produced 
by younger, i.e. more massive progenitors than 
fainter PNe. The h ydrod ynamical simulations by 


ra.mt er _ri\e. me n yaroa ynamicai simulations b} 
Schonberner et al. (120071 see also iMendez et al 


20081 and' Ciardulloll2010l) support such a view and 


provide further c onstraints to the ste llar masses of 
our target PNe. Schonberner et al. 020071 ) found 


that only PNe with central stars masses >0.6 M 0 
can attain the [O III] luminosity of the bright 
cutoff of the PNLF if accompanied by sufficiently 
delayed optically thin/thick transition of the neb¬ 
ular gas. Adopting the empirical i nitial- to-final 
mass relationships of ICatalan et al. ( 2008 ). these 
relatively high core masses imply progenitors with 
an initial mass larger than ~2 M 0 . Our new tar¬ 
gets are between 1 and 2 magnitudes below the 
PNLF cutoff. These slightly fainter magnitudes 
can either correspond to the most luminous PNe 
that have started to fade, or by PNe from slightly 


less massive progenitor s at the their maximum 
luminosity. Figure 7 in Schonberner et a]. (12007 ) 
shows for instance that PNe with core masses of 
0.565 M 0 (initial masses ~1.5 M 0 ) are not able 
to reach the [O III] luminosities of our targets. 
We adopt this value as the lower limit for the PN 
progenitors studied in this work. On the other 
side, it should be considered that with increasing 
mass the number of progenitors decreases (accord¬ 
ing to the standard IMF), and the duration of 
the high [O III] luminos ity phase of t heir nebula e 


becomes much shorter (( Schonber ner et al. 2007). 


Therefore, statistically, it is unlikely to find very 
massive progenitors among our targets. A rough 
upper limit of ~2.5 M 0 may be assumed, consid¬ 
ering that none of them is a type I PN ( Phillips! 
20011 ). We conclude that the PNe presented in 


this work are expected to be produced by stars 
with masses roughly between 1.5 M 0 and 2.5 M 0 . 
The lifetime of these stars is <3.5 Gyr. This mass 
range is also consistent with the modest nitrogen 
enhancement of the target nebulae, becaus e the 


solar-metallicity models of Karakasl ( 201 0) pre¬ 
dict significant N overabundances only for signifi¬ 
cantly more massive progenitors (>4 M 0 ) such 
as o bserved in some G alactic bipolar nebulae 
( Corradi fc Schwarz 1995 ). 

It is important to note that our arguments 
are purely based on single-star evolution consid¬ 
erations. They do not tackle the problem that 
the PNLF cutoff magnitude is the same in all 
galaxies, even in older stell ar populations where 


>2 Mm stars ar e not found (jGiardullo et al.ll20f)2 ; 
ICiardullol 12010 ). This may require alternative 
channels of PN production to populate the bright 
end of the PNLF, such as mergers or mass accre¬ 
tion i n inter acting binaries ( Ciardullo et al 1 120051: 
Soker 20061) . In such case, our mass estimates 


would not be valid. However, low-mass progeni¬ 
tors would not be expected to have the high metal¬ 
licity that we measure. In addition, as discussed 
in the next section, there is independent evidence 
that in the outer regions of M31 2 M 0 stars with 
solar metallicity do exist. 

6. Conclusions 

We have identified two PNe at very large galac- 
tocentric distances with solar oxygen abundances 
and radial velocities consistent with that of the ex- 
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tended disk of M31. Their stellar progenitors, ac¬ 
cording to single-star evolutionary theories, are es¬ 
timated to be in the range 1.5-2.5 M©. These two 
additional objects support our previous conclusion 
(paper II) that the luminous PNe found outside 
the disk of M31 trace a burst of star formation that 
occurred <3 Gyr ago. This interpretation is also 


fully c onsistent with the results of lBernard et al 


( 201 5). who detect trace of such a recent starburst 
in all the HST fields studied, which span a similar 
galactocentric distance range as our PNe. A simi¬ 
lar starburst of solar metallicity stars is also found 
throughout the disk of M31 by the Panchromatic 
Hubble Andromeda Treasury (Ben Williams, pri¬ 
vate c ommunication), as w ell as in the outer disk of 
M33 ( Bernard et al. 2012 ). The most likely expla¬ 
nation for the luminous oxygen-rich PNe in M31’s 
exodisk is that the observed, relatively massive 
stars found in the outer regions of M31 are part 
of the thin disk that has been kicked out during a 
recent encounter with M33, and/or following the 
impact with the unident ified proge nitor of th e Gi ¬ 
ant Stellar Stream (e.g. McConnachie et al • 2009 : 
Bernard et al. 2012 :. 20 lH ). 


Three other PNe in the outer regions of M31 
have O/H abundances ^0.4 dex lower, and show 
significant deviations from the kinematics of M31’s 
exodisk. It is possible that they belong to one 
or another of the halo’s substructures discussed 
in the literature. No bona fide PN belonging 
to the smooth, mostl y m etal-p oor, ancient halo 
described by llbata et ahl ( 20141 ) has been found 
yet. Nor would this be expected since the an¬ 
cient PNe will be far fainter than the younger 
ones. Assuming a total luminosity of t he s mooth 
halo component of ^10 -9 L© (Jlbata et al, 120141) 
and a luminosit y-specific PN number of ~10 - ' 
( Buzzoni. Arnaboldi. fc Corradil 120061) . some 10 2 
halo PNe are expected, but only a tenth of them 
within two magnitudes from the bright cutoff of 
the PNLF, which is the luminosity range that 
we have explored so far. The number could 
be even smaller considering that, in old stel¬ 
lar populations, AGB stars with core masses 
<0.55 M© may not produce PNe, either be¬ 
cause the post-AGB evolution of the core is too 
slow to “light up” the nebulae before they dis¬ 
perse, or because stars escape the AGB phase 
( Buzzoni. Arnaboldi. &; Corradi 20061) . 

Given the complexity of the outer regions of 


M31 and the small number of PNe available, our 
interpretation should be considered as speculative. 
However, it fits well into the modern view of a rich 
interaction and merger history for M31, providing 
independent support to the results obtained using 
other classes of stars, and adding complementary 
information about chemical elements that are best 
studied in ionized nebulae. 

The results of this paper are based on observa¬ 
tions made with (1) the Gran Telescopio Canarias 
(GTC), installed at the Spanish Observatorio del 
Roque de los Mucliachos of the Instituto de As- 
troffsica de Canarias, in the island of La Palma 
and (2) the 3.5 m telescope at Apache Point Ob¬ 
servatory in Sunspot, New Mexico. The Apache 
Point Observatory 3.5-meter telescope is owned 
and operated by the Astrophysical Research Con¬ 
sortium. R.L.M.C. acknowledges funding from 
the Spanish AYA2012-35330 grant. B.B., K.B.K., 
and R.B.C.H. are grateful to our institutions 
and to the NSF for support under grants AST- 
0806490, AST-0808201, AST-0806577, respec¬ 
tively. This research has made use of the USNOFS 
Image and Catalogue Archive operated by the 
United States Naval Observatory, Flagstaff Sta¬ 
tion (http: / /www. nofs. navy, mil/data/fchpix /). 

Facilities: GTC (OSIRIS); APO (DIS); SDSS 
SkyServer 
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A. Flux measurements and computed properties 

Table [2] lists the emission-line measurements of all observed PNe. Column entries are as follows: the first 
column lists the ion and wavelength designation of each line; f(A) gives the value of the reddening function, 
normalized to H/l = 0; and the following lines show the measured flux relative to H/3 = 100. At the bottom 
of each column, for each nebula we list the log of the total observed 11 ,3 flux through the spectrograph slit. 


Note that the latter may differ from the magnitude in Merrett et al. ( 2006lh as not all observations were 
obtained in photometric nights and no attempt was done to account for slit slosses. Table [3] contains the 
logarithmic reddening parameter, c(H/3), and the derived N e and T e of each nebula and the corresponding 
diagnostics. Tables H] [S] display their ionic abundances. They also show the temperatures adopted from 
Table [3] that were used to calculate them. Asterisks denote values that were used in the brightness-weighted 
mean values (wm) in lines below them. Also inclu ded for each elements is the derived ionization correction 
factor ( icf ), calculated as described in iKwitter & Henry (12001 1. Table [7] lists the total abundances. Note 
that the ELSA code determines uncertainties in the dereddened line intensities including the uncertainty in 
c(H/3). Uncertainties in temperature and density diagnostics incorporate the errors in their constituent lines, 
and ionic abundance uncertainties include errors in the relevant line intensities, temperature, and density. 
Final abu ndances inco r porat e errors in the ionic abundances, but not on the icf. For more details, see for 
example Milingo et al. ( 2010f ). Note also that the calculation of a reliable sulfur abundance with the adopted 
icf scheme requires the measurement of S 2+ , that is estimated to be significantly more abundant than S" 1 " in 
all targets. Therefore, in the cases where no S 2+ lines could be measured we do not quote the S/H and S/O 
total abundances. 
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Table 2 

Line identification and observed fluxes 


Ident. 

1(A) 

M2507 

M2538 

M2539 

M2541 

M2543 

M2549 

M2566 

M2988 M31-372 

[O II] A 3727 

0.292 

24.8 

61.5 

20.9 

33.1 

26.8 

27.4 

84.9 

27.1 


Heii+H9 A 3835 

0.262 

4.55 


6.03 







[Nc III] A 3869 

0.252 

44.2 

107 

84.9 

61.9 

93.1 

94.9 

108 

68.2 

48.3 

He I+H8 A 3889 

0.247 

13.8 

8.83: 

16.7 

15.1 

16.0 

17.6 

19.5 

13.2 

18.8 

[Nc III]+He A 3968 

0.225 

28.9 

33.0: 

36.4 

36.8 

42.3 

46.5 

51.5 

34.2 

34.1 

Hel+Hell A 4026 

0.209 



2.22 







[S II] A 4071 

0.196 



1.24 



4.76 

3.56 



Hen A 4100 

0.188 


0.116 a 


0.435 a 


0.308 a 

6.26(-2) a 



H<5 A 4101 

0.188 

24.3 

29.3 a 


25.2 a 


25.0 a 

25.8 a 


22.6 

He ii A 4339 

0.124 


0.215 a 

2.27(-2) a 

0.775 a 

0.105 a 

0.557 a 

0.113 a 



H7 A 4340 

0.124 

44.5 

46.4 a 

46.0 a 

43.2 a 

45.0 a 

44.2 a 

43.5 a 

48.6 

42.2 

[O III] A 4363 

0.118 

6.39 

18.6 

12.9 

13.5 

8.33 

16.0 

8.90 

9.96 

10.3 

He 1A 4472 

0.090 

4.30 

7.63:: 

4.91 

4.21 

5.66 

3.79 

5.32 

5.35 

6.05 

He 11 A 4542 

0.072 




1.62 






N III+O II A 4640 

0.048 




1.14 


2.70 




C III+O II A 4650 

0.045 




0.901 






He 11 A 4686 

0.036 


12.1: 

0.971 

32.5 

4.62 

23.8 

4.95 



He i+[Ar IV] A 4711 0.030 




2.44 


4.39 

3.69 



[Ar IV] A 4740 

0.023 


4.41 

2.46 

1.69 

5.48 

4.53 

3.43 



He 11 A 4859 

0.000 


0.488 a 

4.85(-2) a 

1.62 a 

0.226 a 

1.20 a 

0.242 a 



H/3 A 4861 

0.000 

100 

100 a 

100 a 

100 a 

100 a 

100 a 

100 a 

100 

100 

He 1A 4922 

-0.021 

1.10 


1.38 

1.18 






[O III] A 4959 

-0.030 

232 

479 

378 

324 

418 

415 

427 

311 

324 

[O III] A 5007 

-0.042 

697 

1498 

1140 

978 

1270 

1243 

1302 

948 

985 

[Fe VI] A 5147 

-0.074 



0.421 







Hen A 5412 

-0.134 




2.88 


1.79 




[Cl III] A 5538 

-0.161 



0.466 







[N II] A 5755 

-0.207 

1.24 


0.755 




1.30 



He 1A 5876 

-0.231 

16.8 

18.8:: 

17.8 

12.3 

17.9 

13.7 

17.2 

19.4 

17.0 

[O I] A 6300 

-0.313 

2.41 


2.95 

3.41 




5.26 


Hen A 6311 

-0.315 




0.123 a 






[S III] A 6312 

-0.315 

1.67 


1.08 a 

1.07 a 






[O I] A 6364 

-0.325 

0.753 


1.17 

0.854 




2.15 


[N II] A 6548 

-0.358 

9.52 

33.6:: 

7.35 

8.42 

9.48 

8.48 

27.6 

9.89 

5.22 

Hen A 6560 

-0.360 


1.79 a 

0.146 a 

4.65 a 

0.702 a 

3.68 a 

0.742 a 



Ha A 6563 

-0.360 

326 

33 l a 

311 a 

295 a 

319 a 

318 a 

314 a 

357 

334 

[N II] A 6584 

-0.364 

27.9 

80.0 

21.4 

24.2 

26.6 

26.5 

85.4 

27.8 

15.5 

He 1A 6678 

-0.380 

4.34 


4.60 

3.32 

4.77 

3.92 

4.85 

5.07 

4.79 

[S II] A 6716 

-0.387 

0.758 


0.717 

1.69 

2.22 

1.99 

10.6 

1.03 


[S II] A 6731 

-0.389 

1.56 


1.12 

2.66 

3.00 

2.90 

14.9 

1.50 


He 1A 7065 

-0.443 

10.9 


11.0 

5.29 

7.27 

6.25 

6.13 

11.6 

8.22 

[Ar III] A 7136 

-0.453 

10.4 

17.7:: 

10.8 

5.37 

16.6 

10.0 

22.1 

9.09 

6.38 

He 1A 7281 

-0.475 

1.24 









[O II] A 7324 

-0.481 

18.9 


7.91 

5.58 

4.42 

5.27 

8.61 

11.8 

6.95 

[Ar III] A 7751 

-0.539 

2.77 


2.34 


3.97 

3.66 

4.89 

2.04 


[S III] A 9532 

-0.632 


148:: 








log F(H/3) b 


-15.05 

-15.05 

-15.12 

-15.23 

-15.26 

-15.26 

-15.26 

-15.31 

-15.64 


a Deblended 

b In erg cm -2 s —1 as measured in our extracted spectra 
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Table 3 

Extinction, electron temperatures and densities 



M2507 

M2538 

M2539 

M2541 

M2543 

c(H/3) 

0.17 

0.18 

0.12 

0.06 

0.13 

Te([ 0 III]) 

11300T650 

12700±750 

12250T700 

13150±800 

10250T500 

Te([ N II]) 

16000T4450 

10300 b 

15950T1750 

10300 b 

10300 b 

Te([ 0 II]) 




21950T18100 


Te([ S II]) 






N e ([S II]) 

122f)0 +14750 
±^z,uu_ 12 200 

10000 b 

3000T1600 

3050±1650 

1650±800 


M2549 

M2566 

M2988 

M31-372 


c(H/3) 

0.15 

0.12 

0.29 

0.19 


Te ([0 HI]) 

12950T800 

10350±500 

12250T700 

11950T650 


Ted N HD 

10300 b 

10400±700 

12350±350 a 

12200±350 a 


Te([0 II]) 


14350±5850 




TedS HD 


10250±3500 




M[S HD 

2300T1150 

1900±950 

2250T1100 

10000 b 



a Estimated from T e ([0 III]) according to the prescription from iPagel et al] (jl992i j 
b Default value 
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Table 4 

Ionic Abundances 


Ion 

Tused 


Abundance 


M2507 

M2538 

M2539 

He + 

[0 

III] 

9.85±1.32(-2) 

0.103±0.052 

0.112i0.015 

He+ 2 



1.13±0.36(-2) 

9.04±1.34(-4) 

icf(He) 



1.00 

1.00 

1.00 

0°(6300) 

[N 

II] 

*1.26±0.75(-6) 


*1.48i0.42(-6) 

0°(6363) 

[N 

II] 

*1.23i0.73(-6) 


*1.82i0.53(-6) 

0° 

wm 

1.25i0.74(-6) 


1.57i0.44(-6) 

0+(3727) 

[N 

II] 

*7.56±9.57(-6) 

*4.79±1.01(-5) 

*3.06±1.01(-6) 

0 + (7325) 

[N 

II] 

*1.29i0.78(-5) 


*7.96i3.27(-6) 

o+ 

wm 

9.88±8.30(-6) 


4.41il.25(-6) 

0 +2 (5007) 

[O 

III] 

*1.65±0.43(-4) 

*2.47i0.64(-4) 

*2.08i0.53(-4) 

0 +2 (4959) 

[0 

III] 

*1.59±0.34(-4) 

*2.29i0.48(-4) 

*2.00±0.41(-4) 

0 +2 (4363) 

[O 

III] 

*1.65±0.43(-4) 

*2.47±0.64(-4) 

*2.08i0.53(-4) 

0+ 2 

wm 

1.63±0.40(-4) 

2.43±0.59(-4) 

2.06±0.49(-4) 

icf(O) 



1.00 

1.11±0.07 

l.OliO.OO 

Ar +2 (7135) 

[0 

III] 

*6.20±1.44(-7) 

*8.14±4.42(-7) 

*5.67il.34(-7) 

Ar +2 (7751) 

[0 

III] 

*6.58±1.69(-7) 


*4.98±1.30(-7) 

Ar+ 2 

wm 

6.28±1.44(-7) 


5.55il.30(-7) 

Ar +3 (4740) 




*3.91±0.72(-7) 

*2.72i0.52(-7) 

icf(Ar) 



1.06T0.05 

1.29±0.09 

1.03i0.01 

Cl +2 (5537) 





*3.29i0.67(-8) 

icf(Cl) 





l.OliO.OO 

N+(6584) 

[N 

II] 

*2.60±1.23(-6) 

*1.16i0.20(-5) 

*1.76i0.44(-6) 

N+(6548) 

[N 

II] 

*2.61±1.20(-6) 

*1.44±0.74(-5) 

*1.78i0.39(-6) 

N+(5755) 

[N 

II] 

*2.60±1.23(-6) 


*1.76i0.44(-6) 

N+ 

wm 

2.61±1.21(-6) 

1.24±0.30(-5) 

1.76i0.41(-6) 

icf(N) 



17.5il2.83 

6.72il.48 

48.2il4.76 

Ne +2 (3869) 

[O 

III] 

*3.11i0.77(-5) 

*5.01il.22(-5) 

*4.33il.04(-5) 

Ne +2 (3967) 

[0 

III] 

6.68il.62(-5) 

5.05±1.88(-5) 

6.12il.45(-5) 

icf(Ne) 



1.05i0.06 

1.33±0.11 

1.02i0.01 

S+ 

[N 

II] 

*7.32i7.97(-8) 


*2.73i0.84(-8) 

S + (6716) 

[N 

II] 

7.31i8.04(-8) 


2.70i0.84(-8) 

S + (6731) 

[N 

II] 

7.33i7.93(-8) 


2.75i0.84(-8) 

S+ 

[S 

II] 




S +2 (6312) 

[0 

III] 

*2.17i0.59(-6) 


*l.lli0.31(-6) 

S +2 (9532) 

[O 

III] 


*3.97±2.23(-6) 


icf(S) 



1.57i0.48 

1.24±0.04 

2.37i0.42 
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Table 5 

Ionic Abundances 


Ion 

Tused 


Abundance 


M2541 

M2543 

M2549 

He+ 

[0 III] 

7.63±1.09(-2) 

0.120T0.015 

8.42±1.18(-2) 

He+ 2 


2.98±0.45(-2) 

4.27±0.65(-3) 

2.21±0.33(-2) 

icf(He) 


1.00 

1.00 

1.00 

0°(6300) 

[N II] 

*4.52±0.74(-6) 



O 0 (6363) 

[N II] 

*3.53±0.58(-6) 



o° 

wm 

4.32±0.68(-6) 



0+(3727) 

[N II] 

*1.38±0.37(-5) 

*1.29±0.30(-5) 

*1.13±0.28(-5) 

0+(7325) 

[N II] 

*2.93±0.80(-5) 

*3.37±0.83(-5) 

*2.81±0.75(-5) 

o + 

wm 

1.60±0.26(-5) 

1.58±0.21(-5) 

1.40±0.19(-5) 

0 +2 (5007) 

[0 III] 

*1.45±0.39(-4) 

*4.13±0.99(-4) 

*1.91±0.50(-4) 

0 +2 (4959) 

[0 III] 

*1.40±0.30(-4) 

*3.94±0.76(-4) 

*1.85±0.39(-4) 

0 +2 (4363) 

[O III] 

*1.45±0.39(-4) 

*4.13±0.99(-4) 

*1.91±0.50(-4) 

0+ 2 

wm 

1.44±0.35(-4) 

4.09±0.91(-4) 

1.89±0.46(-4) 

icf(O) 


1.39±0.08 

1.03T0.01 

1.26±0.05 

Ar +2 (7135) 

[0 III] 

*2.59±0.62(-7) 

*1.29±0.29(-6) 

*4.55il.09(-7) 

Ar +2 (7751) 

[0 III] 


*1.24±0.31(-6) 

*6.68±1.77(-7) 

Ar +2 

wm 


1.28±0.28(-6) 

5.12il.22(-7) 

Ar +3 (4740) 


*1.56±0.31(-7) 

*1.03±0.19(-6) 

*4.47±0.86(-7) 

icf(Ar) 


1.50±0.09 

1.07T0.01 

1.34±0.06 

Cl +2 (5537) 





icf(Cl) 





N+(6584) 

[N II] 

*3.45±0.60(-6) 

*4.49±0.71(-6) 

*3.53±0.61(-6) 

N+(6548) 

[N II] 

*3.52±0.61(-6) 

*4.72±0.73(-6) 

*3.33±0.57(-6) 

N+(5755) 

[N II] 




N+ 

wm 

3.47±0.57(-6) 

4.55±0.65(-6) 

3.48±0.57(-6) 

icf(N) 


13.9±2.53 

27.8i5.59 

18.3i3.28 

Ne +2 (3869) 

[O III] 

*2.39±0.59(-5) 

*9.43i2.18(-5) 

*4.04±0.99(-5) 

Ne +2 (3967) 

[0 III] 

4.70±1.14(-5) 

1.41i0.32(-4) 

6.51±1.57(-5) 

icf(Ne) 


1.52±0.11 

1.07i0.01 

1.34±0.07 

S+ 

[N II] 

*1.09±0.28(-7) 

*1.32i0.25(-7) 

*1.05i0.24(-7) 

S+(6716) 

[N II] 

1.10±0.28(-7) 

1.32i0.25(-7) 

1.06±0.24(-7) 

S+(6731) 

[N II] 

1.08±0.28(-7) 

1.32i0.25(-7) 

1.04±0.24(-7) 

S+ 

[S II] 




S +2 (6312) 

[0 III] 

*8.88±2.65(-7) 



S +2 (9532) 

[0 III] 




icf(S) 


1.44±0.08 

1.77i0.17 

1.57i0.11 
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Table 6 

Ionic Abundances 


Ion 

Tused 


Abundance 


M2566 

M2988 

M31-372 

He + 

[O III] 

0.116T0.015 

0.115±0.015 

9.73±1.28(-2) 

He+ 2 


4.57i0.69(-3) 



icf(He) 


1.00 

1.00 

1.00 

0°(6300) 

[N II] 


*4.05±0.82(-6) 


0°(6363) 

[N II] 


*5.13±1.06(-6) 


0° 

wm 


4.36±0.84(-6) 


0+(3727) 

[N II] 

*3.98±1.20(-5) 

*7.50±1.89(-6) 


0 + (7325) 

[N II] 

*6.03±2.59(-5) 

*2.58±0.85(-5) 

*1.16±0.32(-5) 

o+ 

wm 

4.17±1.16(-5) 

1.30±0.27(-5) 


0 +2 (5007) 

[O III] 

*4.10±0.99(-4) 

*1.72±0.44(-4) 

*1.96i0.50(-4) 

0 +2 (4959) 

[O III] 

*3.90±0.75(-4) 

*1.64±0.34(-4) 

*1.88i0.38(-4) 

0 +2 (4363) 

[O III] 

*4.10±0.99(-4) 

*1.72±0.44(-4) 

*1.96±0.50(-4) 

0+ 2 

wm 

4.05±0.90(-4) 

1.70±0.41(-4) 

1.94i0.46(-4) 

icf(O) 


1.04±0.01 

1.00 

1.00 

Ar +2 (7135) 

[O III] 

*1.70±0.38(-6) 

*4.05±0.96(-7) 

*3.34i0.77(-7) 

Ar +2 (7751) 

[O III] 

*1.52±0.38(-6) 

*3.56±0.93(-7) 


Ar+ 2 

wm 

1.67±0.37(-6) 

3.96±0.93(-7) 


Ar +3 (4740) 


*6.22±1.13(-7) 



icf(Ar) 


1.14±0.04 

1.08T0.01 

1.06±0.01 

Cl +2 (5537) 





icf(Cl) 





N+(6584) 

[N II] 

*1.43±0.35(-5) 

*2.65±0.52(-6) 

*1.84i0.35(-6) 

N+(6548) 

[N II] 

*1.36±0.29(-5) 

*2.79±0.54(-6) 

*1.82±0.34(-6) 

N+(5755) 

[N II] 

*1.43±0.35(-5) 



N+ 

wm 

1.41±0.33(-5) 

2.69±0.49(-6) 

1.83±0.33(-6) 

icf(N) 


11.1±3.23 

14.0il.87 

17.8i3.62 

Ne +2 (3869) 

[O III] 

*1.04±0.24(-4) 

*3.86±0.93(-5) 

*2.82±0.68(-5) 

Ne +2 (3967) 

[O III] 

1.64±0.37(-4) 

6.33il.50(-5) 

6.55il.54(-5) 

icf(Ne) 


1.14±0.04 

1.04±0.01 

1.00 

S+ 

[N II] 

*6.58±1.71(-7) 

*3.91i0.94(-8) 


S + (6716) 

[N II] 

6.58±1.72(-7) 

3.91i0.94(-8) 


S + (6731) 

[N II] 

6.58±1.70(-7) 

3.91i0.95(-8) 


S+ 

[S II] 

6.77±5.06(-7) 



S +2 (6312) 

[O III] 




S +2 (9532) 

[O III] 




icf(S) 


1.35±0.08 

1.60±0.11 
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Table 7 

Total Elemental Abundances 


Element 

M2507 

M2538 

M2539 

M2541 

M2543 

He/H 

9.85±1.32(-2) 

0.115T0.053 

0.113T0.015 

0.106T0.012 

0.125T0.015 

N/H 

4.56±1.67(-5) 

8.36±3.26(-5) 

8.50±2.34(-5) 

4.81±1.43(-5) 

1.26±0.37(-4) 

N/O 

0.264T0.107 

0.259T0.089 

0.400T0.056 

0.216T0.047 

0.288T0.057 

O/H 

1.73±0.44(-4) 

3.22±0.71(-4) 

2.12±0.50(-4) 

2.22±0.51(-4) 

4.39±0.95(-4) 

Ne/H 

3.26±0.85(-5) 

6.64±1.63(-5) 

4.43±1.06(-5) 

3.64±0.93(-5) 

1.01±0.23(-4) 

Ne/O 

0.188T0.031 

0.206T0.033 

0.209T0.034 

0.164T0.027 

0.229T0.036 

S/H 

3.52±1.37(-6) 

4.93±2.83(-6) 

2.69±1.05(-6) 

1.44±0.46(-6) 



S/O 

2.04±0.81(-2) 

1.53±0.88(-2) 

1.27±0.33(-2) 

6.46±1.58(-3) 



Cl/H 



3.32±0.67(-8) 




Cl/O 



1.56±0.23(-4) 




Ar/H 

6.66±1.60(-7) 

1.55±0.58(-6) 

8.51±1.61(-7) 

6.22±1.05(-7) 

2.49±0.40(-6) 

Ar/O 

3.85±0.81(-3) 

4.81±1.78(-3) 

4.01±0.63(-3) 

2.80±0.43(-3) 

5.66±0.80(-3) 


M2549 

M2566 

M2988 

M31-372 

Solar ref. 

Orion ref. 

He/H 

0.106T0.013 

0.121T0.015 

0.115T0.015 

9.73±1.28(-2) 

8.50(-2) 

9.70(-2) 

N/H 

6.36±1.89(-5) 

1.57±0.47(-4) 

3.77±0.88(-5) 

3.26±0.70(-5) 

6.76(-5) 

5.37(-5) 

N/O 

0.248T0.048 

0.339T0.077 

0.206T0.032 

0.158T0.026 

0.138 

0.100 

O/H 

2.57±0.59(-4) 

4.64±0.97(-4) 

1.83±0.43(-4) 

2.06±0.48(-4) 

4.89(-4) 

5.37(-4) 

Ne/H 

5.40±1.35(-5) 

1.19±0.27(-4) 

4.03±0.97(-5) 

2.82±0.68(-5) 

8.51(-5) 

1.12(-4) 

Ne/O 

0.210T0.034 

0.255T0.041 

0.220T0.036 

0.137T0.022 

0.174 

0.209 

S/H 





1.32(-5) 

1.66(-5) 

S/O 





2.70(-2) 

3.09(-2) 

Cl/H 





3.16(-7) 

2.88(-7) 

Cl/O 





6.46(-4) 

5.36 (-4) 

Ar/H 

1.28±0.21(-6) 

2.62±0.45(-6) 

4.26±1.00(-7) 

3.54±0.83(-7) 

2.51(-6) 

4.17(-6) 

Ar/O 

4.99±0.70(-3) 

5.64±0.92(-3) 

2.33±0.48(-3) 

1.72±0.36(-3) 

5.13(-3) 

7.77(-3) 
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